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^r 1. Introduction In this paper, we turn our attention to the spin-orbit align- 
2^ ' ment of stars at the bottom end of the initial mass function, where 
Cn Investigation of the spm and orbital properties of solar-type bi- ^^^^^^^ j^ j^^^ common (Lada 2006; Raghavan et al. 2010). In 
T: nanes has been underway for several decades (Weis 1974; Abt ^^e years following the detection of the first brown dwarf Gl 
> & Levy 1976; Bodenheimer 1978; Fekel 1981; Hale 1994, and 229B by Nakajima et al. (1995), a number of surveys yielded 
• --H references therein). A large fraction of close ( <40 AU) solar- ^^e discovery of a large number of low mass star ( < 0. 1 M^) and 
^ type binaries have been found to have coplanar spms and or- ^rown dwarf binary systems, e.g. Burgasseret al. (2007, and ref- 
^ Dits (Hale 1994). Studies by Weis (1974) yielded similar re- grences therein). Following these discoveries, the introduction of 
- - -suits over a wider range of spectral classes [B - F]. More recent j^^^j. g^j^g ^^^^. (lgS) adaptive optics (AO) systems on ground- 
studies have turned up a diversity of configurations. Jensen et ^ased telescopes provided the means of assessing the dynami- 
al. (2004) and Monm et al. (2006) have reported planar align- ^^j ^^^^ ^f ^^^^ ^^^^^^^ ^^^^^ ^^ ^j 2004; Dupuy et al. 2010; 
ment and misalignment for wider-sepai-ation binaries. Similar Konopacky et al. 2010). More recently, Konopacky et al. (2012) 
to these studies, in the case of very close binaries of semi-major obtained resolved LGS AO spectroscopic measurements of indi- 
axis -0.3 AU, there have been some examples where systems ^j^^^j component rotation velocities for a sample of eleven very 
exhibit both ahgned and misaligned axes (Albrecht et al. 2009, jo^ mass (VLM) dwarf binai'ies (M,„, < 0.185 M^, as defined by 
201 1). More recently. Wheelwright et al. (201 1) have also re- ^lose et al. (2003)). These data provided additional parameters 
ported coplananty between HAe/Be binaiy systems and circum- f^^. ^^^^^^^ separated by ~1 - 10 AU, but could only be used to 
stellar disks. Thus, it appears that while the spin-orbit alignment tentatively investigate spin-orbit alignment, since other parame- 
is common, there are exceptions at all orbital separations. ^^^^ s^^h as individual component rotation periods, and system 



ABSTRACT 

Studies of solar-type binaries have found coplanarity between the equatorial and orbital planes of systems with <40 AU separation. 
By comparison, the alignment of the equatorial and orbital axes in the substellar regime, and the associated implications for formation 
theory, are relatively poorly constrained. Here we present the discovery of the rotation period of 3.32 ± 0.15 hours from 2MASS 
J0746+20A - the primary component of a tight (2.7 AU) ultracool dwarf binary system (L0+L1.5). The newly discovered period, 
together with the established period via radio observations of the other component, and the well constrained orbital parameters and 
rotational velocity measurements, allow us to infer alignment of the equatorial planes of both components with the orbital plane of the 
system to within 10 degrees. This result suggests that solar-type binary formation mechanisms may extend down into the brown dwarf 
mass range, and we consider a number of formation theories that may be applicable in this case. This is the first such observational 
result in the very low mass binary regime. In addition, the detected period of 3.32 ± 0. 15 hours implies that the reported radio period 
of 2.07 + 0.002 hours is associated with the secondary star, not the primary, as was previously claimed. This in turn refutes the 
claimed radius of 0.78 ±0.lRj for 2MASS J0746+20A, which we demonstrate to be 0.99 + 0.03 Rj. 

Key words, binaries: close - brown dwarfs - stars: formation - stars: low-mass - stars: rotation 
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properties inferred from evolutionary models, such as radii, were 

led. 

Article number, page 1 of 9 



Table 1. Properties of 2MASSW J0746425+200032AB 



Parameter 



Rotation period (hrs) 

V sin / (km s"') 

Equatorial Velocity (km s"') 

Period ratio 

V sin i ratio 

Orbital period (yrs) 

Semi-major axis (mas) 

Inc. ORB (deg)* 

Inc.EQ (deg)^ 



2MASS J0746+20A 
^015^ 



2MASS J0746+20B 



3.32: 

19 + 2 

36 + 4"l" 
62 +""2 

^■-^'-0.10 

12.71 + 0.07 

237 3+'-^ 
41.8 ±0.5 

32 + 4"'' 



2.07 + 
33 

56: 



0.002 

t3 

:2t 



12.71 ±0.07 

237 3+^-^ 

41.8 ±0.5 

36 ±4^ 



Age (log yrs) 


9.1 ±0.1' 


9.1 ± 0.1 ' 


MasStoto, (Mo) 


0.151 ±0.003 


0.151 ±0.003 


Mass (Mq) 


0.078 ± 0.004"'" 


0.073 ± 0.004"'" 


Lithium? 


No 


No 


Radius (Rj) 


0.99 ± 0.03^ 


0.96 ± 0.02' 


Gravity (log g) 


5.34 ± 0.02^ 


5.34 ± 0.02"'" 


hhoi (log LILq) 


-3.64 ± 0.02 


-3.77 ± 0.02 


Abs. mag (7) 


11.85 ±0.04 


12.36 ±0.10 


Abs. mag (//) 


11.13 ±0.02 


11.54 ±0.03 


Abs. mag (K) 


10.62 ± 0.02 


10.98 ± 0.02 


References 


1,3-5 


1-5 



Notes. ^ Rotation period discovered in this work; all other parameters derived here based on evolutionary models of Chabrier et al. (2000). 
Xlnc.oRB is the system orbital inclination as measured by Konopacky et al. (2012), whereas lnc.£^ is the equatorial inclination of each component, 
calculated in this work. 
References. (1) This work; (2) Berger et al. (2009); (3) Konopacky et al. (2010); (4) Bouy et al. (2004); (5) Chabrier et al. (2000). 



Previous studies have shown that the presence of mag- 
netic fields can affect binary formation (Mestel 1977; Boden- 
heimer 1978; Fekel 1981; Li et al. 2004, and references therein), 
whereby such fields could potentially contribute to a loss of an- 
gular momentum on large scales, or the tilting of stellar spins on 
small scales. Are these concerns of special relevance in the ultra- 
cool dwarf regime? M dwarfs later than M3 are now associated 
with intense magnetic activity, often possessing surface mag- 
netic field strengths of a few kG and greater (Reiners & Basri 
2007; Hallinan et al. 2008; Berger et al. 2009; West et al. 201 1). 
In this paper, we now have sufficient data to investigate the or- 
bital properties for the magnetically-active VLM dwarf binary, 
2MASS J0746+20AB (M,,,, = 0.151 ± 0.003 Mq). This detec- 
tion can shed light on the formation of VLM binary stars, and 
could signal that a scaled-down version of the formation mecha- 
nism for solar-type binary systems holds in this regime, despite 
the presence of a ~1.7 kG magnetic field in this case (Antonova 
et al. 2008). Characterizing the fundamental properties of VLM 
binary star formation is important in establishing a correlation, 
if any, in the formation and evolution of all types of binary stars. 



2. 2MASSW J0746425+200032AB 

2.1. Properties 

2MASS J0746+20AB is an L dwarf binary (L0+L1.5) that is 
located at a distance of 12.20 ± 0.05 pc (Dahn et al. 2002). It 
is a tight binary system, with a separation of ~2.7 AU (Reid 
et al. 2001) and an effective temperature of between 1900 - 
2225 K (Vrba et al. 2004). The latest high-precision dynamical 
mass measurements yielded a total mass of 0.151 ± 0.003 M© 
(Konopacky et al. 2010), initially measured to be 0.146 ±ooog 
Mq by Bouy et al. (2004), placing the dwarf in the VLM bi- 
nary regime. Bouy et al. (2004) based their individual compo- 



nent masses on a model-derived age estimate. However, Gizis 
& Reid (2006) questioned these individual mass estimates, and 
argued against the secondary component being a brown dwarf as 
predicted by Bouy et al. (2004), and instead favored a substellar 
or low mass star classification. 

The L dwarf was reported as an active radio source by 
Antonova et al. (2008), who estimated magnetic field strengths 
of ~ 1.7 kG based on the detection of a single highly polarized 
burst of emission. The dominant emission was quiescent, the na- 
ture of which is still debated. Following this observation, Berger 
et al. (2009) reported radio emission with a rotation period of 
2.07 ± 0.002 hours. Simultaneously, they detected periodic Ha 
emission. The period was the same in both instances, and was 
consistent with stellar rotation. They also estimated a stellar ra- 
dius of 0.078 ± 0.010 Ry for 2MASS J0746+20A; in order to 
do this, they assigned an average v sin / of 27 ± 3 km s"' (taken 
from a number of previous studies (Reid et al. 2002; Bailer-Jones 
2004; Reiners & Basri 2008)) to the primary. This assumption 
was based on the observations of Bouy et al. (2004), who estab- 
lished that 70% of the light contribution was coming from the 
primary star. Thus, they argued that both the Ha and radio emis- 
sion were also emanating from the primary component of the 
system. In addition, marginal evidence of periodicity (of a few 
hours) has previously been reported by Clarke et al. (2002) and 
Bailer-Jones (2004). Target properties are shown in Table 1 . 



2.2. Instrumentation and observations 

We obtained observations over a ~2 year baseline for 2MASS 
J0746+20 to search for periodic photometric variability from ei- 
ther, or both components of the binary system. The observa- 
tions encompassed 4 separate epochs in January 2009, February 
2010, November 2010 and December 2010. These were carried 
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Table 2. Observation details 



Source 


Date 


Length 


Exp. 


Photometric 


Band 


Telescope 




of obs. 


of obs. 


time 


error 




/ instrument 




(UT) 


(~hrs) 


(s X CO add) 


(%) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


2MASS J0746-t-20AB 


2009 Jan 25 


6.0 


25x1 


0.21 




VATT/4K 




2009 Jan 26 


6.8 


25x1 


0.28 




VATT/4K 




2009 Jan 28 


7.4 


25x1 


0.24 




VATT/4K 




2010 Feb 19 


4.5 


5x 12 


0.27 




VATT/GUFI 




2010 Feb 20 


4.0 


5x 12 


0.30 




VATT/GUFI 




2010 Nov 13 


4.6 


5x 12 


0.31 




VATT/GUFI 




2010 Nov 14 


5.5 


5x12 


0.33 




VATT/GUFI 




2010 Dec 2 


6.0 


5x12 


0.25 




VATT/GUFI 




2010 Dec 12 


3.0 


5x12 


0.29 




VATT/GUFI 




2010 Dec 13 


6.8 


5x 12 


0.32 




VATT/GUFI 




2010 Dec 14 


7.0 


5x 12 


0.34 




VATT/GUFI 



Notes. Column (1) Target of campaign. (2) Date of each observation in UT. (3) The total time observed per night for each observation, in hours. 
(4) The exposure time for each observation in seconds, followed by the binning factor used to increase data point S/N. (5) The mean photometric 
error for light curve data points per night. (6) The VATT Johnson photometric waveband used (I: ~7200 - 9100 A). (7) Telescope and detector 
used for a given observation. These observations were taken over 4 separate epochs, spanning ~62 hours of data over a 2 year baseline. 



out with the GUFI (Galway Ultra-Fast Imager) and VATT 4K 
photometers on the VATT telescope^ using the VATT I-Arizona 
filter (~7200 - 9100 A). GUFI was first commissioned by as- 
tronomers at NUI Galway to operate on the 1.5 m Loiano tele- 
scope in Bologna, Italy (Sheehan & Butler 2008). It was mod- 
ified thereafter for use on the VATT telescope, where it is cur- 
rently stationed as a visitor instrument. GUFI was developed as 
a dedicated high-speed photometer for a long-term observational 
campaign that focuses on the periodic variability of a sample of 
radio detected ultracool dwarfs, and the associated mechanisms 
responsible for such periodic signals (Harding et al. 2013). It 
provides a FOV of ~ 3' x 3', a plate scale of 0.35" pixel"^ and 
performs high-time resolution imaging (e.g. 34 frames per sec- 
ond full frame, with ~2 ms readout rates). Data was taken with 
exposure times of 5 seconds where frames were later summed 
in image space to 1 minute to increase the signal to noise (S/N). 
GUFI was used for the February 2010, November 2010 and De- 
cember 2010 epochs. 

We used the VATT 4K photometer for the first epoch of the 
2MASS J0746+20 observations in January, 2009, before GUFI's 
commissioning in May 2009. The VATT 4K is the facility pho- 
tometer of the VATT observatory, and has a native plate scale of 
0.188" pixel-' and a FOV of -12.5' x 12.5'. Details of obser- 
vations, and the typical photometric error per given observation 
are shown in Table 2, where the full sample includes ~62 hours 
of observations. 



pending). We also carried out photometry for all reference stars 
in order to measure their level of variability - this comparison 
ensured that the periodicity was indeed intrinsic to the binary. 
Reference stars were chosen after assessing their stability (in 
terms of flux) throughout the night, their isolation on the CCD 
chip, the properties of their seeing profiles, and their color in- 
dex with respect to the target. Photometric apertures (in pixels) 
which provided the highest S/N for the target star were selected 
for aperture photometry. Since GUFI did not resolve each com- 
ponent of the binary system, the chosen aperture included the 
combined flux from both stars. Differential photometry was ob- 
tained by dividing the target flux by the mean flux of selected 
reference stars. Finally, after differential photometry was carried 
out, we assessed the photometric light curves for evidence of 
any residual systematic trends due to e.g. the effects of increas- 
ing/decreasing airmass. These effects were negligible through- 
out the campaign, where typical seeing at the VATT observatory 
was ~0.7 - 1.6". However, where small trends were observed 
due to such second order extinction effects, we fitted a linear (or 
quadratic, trend depending) fit to the data and then applied an 
inverse form of the chosen fit. Removing trends is important in 
a variable light curve, since they can introduce aliasing effects 
in the subsequent spectral analysis (discussed in Section 3). The 
calculated photometric errors are shown in Table 2 (column 5), 
and the procedure for calculating these errors is outlined in Sec- 
tion 3. 



2.3. Data reduction and differential photometry 

Standard data reduction techniques were employed via a GUFI 
pipeline (Sheehan & Butler 2008), which performs bias subtrac- 
tion (using zero-integration frames) and flat-fielding (using twi- 
light flat-fields). The FOVs of GUFI and the VATT 4K provided 
up to 20 reference stars for the target field (detector FOV de- 



' The Vatican Advanced Technology Telescope (VATT) tele- 
scope (1.83 m) facility is operated by the Vatican Obser- 
vatory, and is part of the Mount Graham International Ob- 
servatory. Details of VATT detectors can be found here: 
http://vaticanobservatory.org/VATT/index.php/telescope-instruments/. 



3. Periodic Variability 

3. 1 . Variability analysis 

The Lomb-Scargle (LS) periodogram (Lomb 1976; Scargle 
1982) was used to search for periodic variability in the binary 
light curves. This technique is very effective for unevenly spaced 
data, and uses the discrete Fourier Transform (DFT). The algo- 
rithms output power spectra with a range of peaks of different 
amplitudes. The significance of these peaks are then analyzed to 
a given sigma level as computed by the LS algorithm, which cor- 
responds to possible periodic variability. We selected a range of 
peaks of >5cr significance, and then inspected the various solu- 
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HJD: 2454856.71001 1647 



HJD: 2454857.631957755 



HJD: 2454859.653165393 





January 25 2009 UT (hrs) 



4 6 8 10 

January 26 2009 UT (hrs) 




4 6 8 10 

January 28 2009 UT (hrs) 



HJD: 2455246.705281595 



HJD: 2455247.685167070 



HJD: 2455513.820684682 




5 6 7 8 9 

February 19 201 OUT (hrs) 




5 6 7 

February 20 2010 UT (hrs) 




9 10 11 12 

November 13 2010 UT (hrs) 



HJD: 2455514.814986587 



HJD: 2455532.794271032 



HJD: 2455542.919808880 




9 10 11 12 

November 14 2010 UT (hrs) 




10 12 

December 2 2010 UT (hrs) 




11 12 

December 12 2010 UT (hrs) 



HJD: 2455543.815873595 



HJD: 2455544.785973326 



Reference star lightcurve 




9 10 11 12 13 

December 13 2010 UT (hrs) 




10 12 

December 14 2010 UT (hrs) 




9 10 11 12 

UT (hrs) 



Fig. 1. I-band light curves of 2MASS J0746+20AB. Each figure shows relative flux vs. UT times (and dates). We also mark the HJD time above 
each figure corresponding to the start-point of each observation. All of the data presented here was taken in the ' VATT Arizona' I-band broadband 
filter (~7200 - 9100 A), where the baseline extends over ~2 years - as shown in Table 1. In this work, we report the discovery of consistent periodic 
variability from 2MASS J0746+20A, with a period of 3.32 ± 0.15 hours, and an amplitude variability range of 0.40 - 1.52%. We note that the 
full in-depth analysis investigating the stability of the target's phase and amplitude is discussed in other work (Harding et al. 2013). Some poor 
weather conditions during constituent epochs were present - e.g. January 25 & 26 2009 were taken under thin cloud and high wind conditions 
and thus were binned by a factor of 2 compared to other data sets (Table 1). The vertical arrow marked on the November 14 2010 light curve 
denotes a period of complete cloud cover, and so data was removed accordingly. The red light curve (bottom right) is an example of a reference 
star plotted against all others (taken from November 13 2010), illustrating the stability of the chosen reference stars with respect to the target. The 
mean amplitude variability of the reference stars is shown in Table 3. 



tions by phase folding the Hght curves to these periods. Compar- 
ing these phase folded light curves from both individual nights, 
and from different epochs, allowed us to assess the level of agree- 
ment of phase, for a given period. 

The standard deviation of the phase fold was also considered. 
We calculated this via Phase Dispersion Minimization (PDM) 
techniques (Stellingwerf 1978). PDM is a least squares fit cal- 
culation where the coiTect period used in the routines produces 
the least data point scatter in the resulting phase fold. This cal- 



culation is called the 'PDM Theta statistic' (0), where the min- 
imum value is the most likely period solution. We selected a 
broad range of periods (e.g. 1-5 hours), and ran the PDM al- 
gorithm giving a corresponding range of minima. The signifi- 
cance of each was assessed via 10^ Monte Carlo trials, which 
tests whether any detected minimum could be a result of noise 
alone. 

Finally, we assessed the peak to peak (PtP) amplitude vari- 
ability of the target light curves (PtPrargei) by fitting a model sinu- 
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Table 3. Periodic photometric variability from 2MASS J0746+20A 



Phase fold of all epochs 



Parameter 


2M J0746+20A 


2M J0746-H20B 


Rotation period (hrs) 


3.32 + 0.15 


2.07 + 0.002 


LS ' period (hrs) 


3.318 




LS sign.l (a-) 


»5 




PDM* period (hrs) 


3.32 




PtPmrge, (%) 


0.40-1.52 




Mean (Treference (%) 


-0.36 




Variability 


Photometric 


Radio/Ha 


References 


1 


2 



Notes. ^Period obtained from the Lomb-Scargle (LS) periodogram anal- 
ysis, as shown in Figure 3 (left). We also include the significance/false 
alarm probability of the detection as calculated by the LS algorithm. 
* Period obtained from the Phase Dispersion Minimization (PDM) tech- 
nique, shown in Figure 3 (right). 

References. (1) This work; (2) Berger et al. (2009). 



soidal signal to the raw data using the detected period, and varied 
the phase and amplitude of the model. A Chi squared (x^) min- 
imization was performed giving the best fit for a given night's 
peak to peak amplitude variability. The scatter in the reference 
star light curves ((Treference) Were established by calculating their 
standard deviation. 



3.2. Photometric error and period uncertainty estimation 

We estimated the error in the relative magnitude of the target star 
using the iraf.phot^ task. The error in magnitude was then con- 
verted to an error in flux (since we plot relative flux vs. UT in 
our photometric light curves). By using iraf.phot, we take both 
formal (e.g. flat-fielding) and informal (e.g. fringing^) errors in 
to account - these are usually quite difficult to assess indepen- 
dently. 

The period uncertainty was calculated for individual nights 
using a.x^ test. We could not run this test on the entire baseline 
at once due to large gaps in the time series. Furthermore, we 
did not achieve an accurate enough period in consecutive epochs 
to phase connect the 2 year baseline together The quoted pe- 
riod and eiTor in this work was therefore derived within epoc hs, 
which was found to be consistent for all. 



3.3. Results 

Although we do not resolve each component of the binary as a 
point source, most intriguingly, we show optical periodic modu- 
lation of 3.32 + 0. 15 hours in VATT I-band (Figure 1). We show 
phase folded light curves for this period in Figure 2. Therefore, 
this optical periodic variability originates from the other compo- 
nent to that producing the radio emission - reported by Berger 
et al. (2009) where the binary exhibited periodic bursts of ra- 
dio emission of 2.07 + 0.002 hours. By adopting the estimated 
radii in this work (discussed in the next section), in addition to 
the V sin / measurements shown in Table 1 (Konopacky et al. 
2012), we derive maximum period values of ~4.22 hours and 



Image Reduction and Analysis. Found here: http://iraf.noao.edu/. 
' Fringing is an additive optical effect at red/NIR wavelengths that oc- 
curs in the thinned substrate of back-illuminated CCDs. This is due 
to atmospheric spectral emission such as OH, and varies in amplitude 
across the frame. Thus it needs to be removed if varying at a level 
greater than the amplitude variability of the target star. In the case of 
2MASS J0746-I-20, we assessed this effect and found it to be negligible. 



Raw 

Bin factor 5 



.•-.'< 



.•-.'* 










Fig. 2. Phase folded light curves of all epochs from the 2MASS 
J0746+20A observation campaign. The black light curve is folded us- 
ing the raw light curves as shown in Figure 1, and the red is the same 
binned to a factor of 5. Each light curve is phase folded to the detected 
period of 3.32 ± 0.15 hours. 



-2.38 hours for 2MASS J0746-I-20A and 2MASS J0746H-20B, 
respectively. This infers that we detect the optical rotation pe- 
riod of 2MASS J0746-I-20A, the primary component of the sys- 
tem. This result is contrary to what was claimed by Berger et 
al. (2009), who attributed the radio period to the primary based 
upon an incoiTect assumption that the v sin / broadening was due 
to the primary star, since it is more luminous (Bouy et al. 2004). 

A LS periodogram showing the significance of the detection 
in our work is shown in Figure 3 (left), and the PDM analy- 
sis in Figure 3 (right). Both analyses support a period of 3.32 
hours for one binary component. The PtPiargei variability range 
for all light curves was found to be ~0.40 - 1.52%, where we 
also observe changes in light curve morphology throughout the 
campaign. The source of this behavior is discussed in other work 
(Harding et al. 2013). However, in this paper we do investigate 
the presence of the other binary component - since variable chro- 
mospheric emission and an active magnetosphere was reported 
(Berger et al. 2009). Thus the presence of magnetic surface fea- 
tures such as hot or cool spots must be considered for both com- 
ponents. We highlight a cluster of reasonably low false alarm 
probability power spectra around ~12 days"' (or ~2 hours), ap- 
proximating the radio period. To fuither assess the possible pres- 
ence of an underlying period from the secondary component, we 
fitted a sinusoidal signal with a period of 3.32 hours to each of 
the raw light curves. We then subtracted this out, and re-ran 
LS periodograms on each of the remaining data points. For each 
night, those probabilities were no longer present and thus did not 
support the presence of the 2.07 + 0.002 hour period. We also 
point out multiple power spectra and minima around the re- 
ported rotation periods in Figure 3 (left and right). These peaks 
and minima are present as a result of spectral leakage, which is 
due to large gaps in the data between consecutive observations. 

Binary system rotation periods are an important parameter in 
the diagnosis of the system's orbital coplanarity. A radius esti- 
mate, in addition to a v sin / measurement of a system can only 
loosely infer spin-orbital alignments - the rotation period is an 
essential variable in the calculation of 'v', which makes its dis- 
covery all the more pertinent. Furthermore, the period reported 
here has allowed us to infer which period matches each binary 
member, and thus allowed for an effective estimate of masses 
and radii, which we discuss in the following section. 
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Power spectrum for 2MASS J0746+2000A (all epochs) 



Phase Dispersion Minimization 




10 15 

Frequency (days~^) 



3.4 3.6 

Period (hours) 



Fig. 3. [LEFT] Lomb-Scargle periodogram of all 2MASS J0746+20A epochs, calculated from the combined dataset from Figure 1. The 
three red dashed-dotted horizontal lines (top-to-bottom) represent a Scr, 3cr and Icr false-alarm probability of the peaks, as determined by the LS 
periodogram algorithm. The x-axis is plotted in days" ' , since each epoch was time-stamped in units of Heliocentric Julian Days (HJD) - for period 
accuracy over such a lengthy baseline. [RIGHT] Phase Dispersion Minimization (PDM) plot of the entire ~2 year observation baseline. Here we 
plot of the PDM 'Theta (©) statistic' vs. Period (in hours), where the routines minimize the value at 3.32 hours. 



4. An estimate of age, mass and radius 

By adopting the established total system mass of 0.151 + 0.003 
Mq, the photometric J H K measurements and bolometric lu- 
minosity measurements of Konopacky et al. (2010), in addition 
to the lack of detected lithium in the binary dwarf's spectrum 
(Bouy et al. 2004), we were able to place constraints on the evo- 
lutionary models of Chabrier et al. (2000), in determining the 
mass range and radii of each component. We make no initial as- 
sumptions for the age of the system, however young ages were 
ruled out based on this absence of lithium. Thus we had three 
measured quantities to estimate the mass track (i.e. J H K col- 
ors, Lhoi, and Li=0). Once we identified a range of ages that did 
not contain lithium for either component, we interpolated over a 
range of masses based on the correlation between the J H K col- 
ors of Konopacky et al. (2010), and those of the Chabrier et al. 
(2000) models, and then over the bolometric luminosities - thus 
establishing the best agreement between each quantity. Most 
importantly, assuming the stars are coeval, the masses were ul- 
timately constrained since the sum of each component's mass 
could not be more than the estimated total mass of 0. 15 1 ± 0.003 
Mq for the system. 

We find an age of ~ 1 - 1.5 Gyr for the binary based on this as- 
sessment, as well as individual mass estimates of 0.078 ± 0.004 
Mo and 0.073 + 0.004 Mq for 2MASS J0746-)-20A and 2MASS 
J0746-H20B, respectively. These mass estimates are largely con- 
sistent with Bouy et al. (2004), and in good agreement with Gizis 
& Reid (2006) and Konopacky et al. (2010), and furthermore 
infer that each component lies at, or just below, the substellar 
boundary, supporting the prediction of a low mass star classifica- 
tion for the secondary member (Gizis & Reid 2006). The differ- 
ence in rotational velocity between these stars is most intriguing. 



considering the similarity of component mass estimates. Per- 
haps there is a larger difference in mass between each member 
than what has been infeiTed from the evolutionary models in the 
above studies. 

By contrast, an age of ~1 - 1.5 Gyr identifies the system as 
a much older binary dwarf than originally predicted by Bouy 
et al. (2004), who found the system to be -150 - 500 Myr old. 
This is a large discrepancy. However, we point out that Bouy et 
al. (2004) put forward these ages despite the absence of lithium 
in the binary's spectrum - which is expected to be present for 
stars of this age. We note however, that some studies suggest 
that lithium can indeed be depleted at younger low mass star 
ages, due to the effect of episodic accretion (Baraffe & Chabrier 
2010). Although the absence/presence of lithium was originally 
used as a test for sub-stellarity (e.g. Rebolo et al. (1992)), it is 
consequently not necessarily a robust indicator of stellar age. A 
much younger age is inconsistent with the surface gravity esti- 
mates of Schweitzer et al. (2001), who compared spectra to the 
models of Allard et al. (2001). Their temperature estimates do 
agree with those of Bouy et al. (2004), but their estimated grav- 
ity is too high for a ~150 - 500 Myr old object. These gravity 
estimates (from high-resolution spectra: log g ~5.0 [K I /17685] 
- 5.5 [Ca I /15673]; and from low-resolution spectra: log g ~ 6.0) 
are more consistent with the ages we infer in this work. 

Therefore, adopting the estimates of age and mass above 
places each star just below 1 Rj, as shown in Figure 4, where 
we estimate radii of 0.99 + 0.03 Rj for 2MASS J0746-H20A and 
0.96 + 0.02 Rj for 2MASS J0746-H20B. These predictions are 
consistent with those of Konopacky et al. (2010). However, the 
radius estimate for 2MASS J0746+20A is inconsistent with the 
estimate of Berger et al. (2009), due to their assumption of a v 
sin / of 27 km s"' for 2MASS J0746+20A, and a period of 2.07 
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Fig. 4. Isochrones of radii (Rj) vs. mass (Mq) for 2MASS 

J0746-I-20AB, for ages of ~1 (log 9.0; red), -1.25 (log 9.1; blue) and 
~1.5 Gyr (log 9.2; green), derived in this work by using the evolution- 
ary models of Chabrier et al. (2000). The errors associated with the 
mass and radii estimates are based on the constraints of the Chabrier 
et al. (2000) models for a given age, as outlined in Section 4. We also 
include evolutionary tracks for ages of ~150 - 500 Myr (log 8.2 - 8.7 
years) as predicted by Bouy et al. (2004). These are shown by the grey 
dash-dotted lines. Based on the work this paper, the binary appears to 
be a much older age than the range predicted by Bouy et al. (2004), 
which require much larger radii for both components. 



+ 0.002 hours. This rotational velocity is approximately 30% 
larger than the established v sin / of 19 + 2 km s"' (Konopacky 
et al. 2012), thus incorrectly placing their radius estimate much 
lower than those predicted by evolutionary models, at 0.78 + 0. 1 
Rj. Our newly established radii are much more consistent with 
model predictions, based on our identifying the primary compo- 
nent as the non-radio pulsing star, and strong evidence support- 
ing a spin-orbit alignment to within 10 degrees for both stars. 
See Table 1 for a summary of system properties, including our 
estimates of surface gravity. 

5. Inferred spin-orbit alignment 

The discovery of the rotation period of 2MASS J0746H-20A 
in this paper, in addition to the spectroscopic observations of 
Konopacky et al. (2012) and the radio observations of Berger 
et al. (2009), yield a v sin / ratio of 0.57^q[q, and a period ra- 
tio of 0.62^QQ3. By adopting these v sin / and rotation period 
measurements, the estimated radii in this work indicate that both 
components in the binary have their equatorial spin axis aligned 
with the orbital plane to within 10° - consistent with the trends 
for solar-type binaries with separations < 40 AU (Hale 1994). 
Is this evidence for a common formation pathway across a wide 
swath of stellar masses? Our best fit indicates that the axes are 
orientated with respect to the observer at 32 + 4° and 36 + 4°, 
respectively (Figure 5 and Figure 6). If the measurements are 
taken at face value, both objects could also be perfectly coplanar, 
implying that evolutionary models over-predict their radii for an 
assumed mass and age. The inclination angle of the equatorial 
spin axis of each component, sin i^ and sin is, were derived as 
follows: 



V sin lA 
2-7T- Fa/Pa 



and sin ig 



V sin iB 
1-n- relPB 



(1) 



where v is the rotational velocity, r is the radius and P is the 
period. We estimated the uncertainly in equatorial inclinations 
above by using Monte Carlo simulations. Each of the variables 
in Equation 1 had known errors as calculated by Konopacky et al. 
(2012) (v sin /), this work (r^ and rg, P^) and Berger et al. (2009) 
(Pb), respectively. For each of these variables, we generated 1 0^ 
copies, where each iteration had random noise (normally dis- 
tributed) added based on the reported errors. For each of the 10^ 
simulations, we measure the standard deviation of the variables 
and then calculated a mean value for each. This mean standard 
deviation over all of the simulations is thus taken as the uncer- 
tainty (randomly calculated) in each quantity and applied to the 
final calculation of the equatorial inclinations in Equation 1 (as 
shown in Table 2 and Figure 5 & Figure 6). 

Finally, we highlight a geometric effect with respect to sys- 
tem spin axes inclinations. The inclinations of the rotation axes 
of each component to our line of sight may be equal, but this 
does not necessarily imply that the orbital planes are coplanar 
For example, two edge-on spins can still be orthogonal on the 
sky, and could be coincidentally equal. 

6. Discussion: implications for formation 

This is the first study to assess the orbital alignment properties 
of a VLM binary system. The observed orbit-spin alignment 
is consistent with several different formation pathways. There 
are numerous binary (and multiple) star system formation the- 
ories, the most prominent of which are turbulent core fragmen- 
tation, disk fragmentation, and formation via dynamical interac- 
tions both during the main accretion phase and after, e.g. (Hoyle 
1953; Adams et al. 1989; LaugMin & Bodenheimer 1994; Bon- 
nell 1994a; Tsuribe & Inutsuka 1999; Padoan & Nordlund 2002; 
Bate et al. 2003; Bonnell et al. 2004; Clark et al. 2008; Kratter 
et al. 2010; Kratter 2011). A nearly coplanar system might re- 
sult from any of these models, although formation via a dynam- 
ical interaction at late times is unlikely; tidal interactions are too 
weak at these distances, so alignment would be coincidental. 

Disk fragmentation naturally produces aligned systems, and 
can also drive components towards equal mass (Adams et al. 
1989; Laughlin & Bodenheimer 1994; Bonnell 1994a; Bonnell 
& Bate 1994b; Kratter et al. 2010). However, disk fragmenta- 
tion is more likely for higher mass systems (>1 Mq) (Kratter et 
al. 2008; Offner et al. 2010). Notably, these models have not 
been extended down to the brown dwarf regime. Stamatellos 
& Whitworth (2009) proposed an alternative disk fragmentation 
scenario, where binary brown dwarfs are born within the disk 
of a more massive star. This scenario is somewhat inconsistent 
with 2MASS J0746+20AB, as it produces only very high eccen- 
tricity binaries (they predict e > 0.6 compared to this system's 
e - 0.487 (Konopacky et al. 2010)). It is also unclear whether 
ejection would preserve spin alignment. 

Both core fragmentation and competitive accretion can also 
produce aligned systems. In the former scenario, fragments may 
share the core's net angular momentum vector (Matsumoto & 
Hanawa 2003). We note that while equal inclinations is a neces- 
sary condition for aligned angular momenta, it is not a sufficient 
one, unless the inclinations are zero. Recent work by Jumper 
& Fisher (2012) has shown that a straightforward extrapolation 
of the turbulent core fragmentation model to lower masses nat- 
urally reproduces the separation distribution of brown dwarfs. 
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Fig. 5. Equatorial rotational velocity (y-axis, left) vs. inclination angle 
of the orbital plane (x-axis, bottom). The measured v sin i of 2MASS 
J0746+20A is shown by the red solid curve. The dashed lines in all 
cases represent the associated errors. The green vertical solid line (x- 
axis, top) highlight the alignment of the spin-orbit axes. Y-axis, right 
corresponds to the radius of the dwarf (Rj = ~69550 km) as calculated 
in this work, where we have marked the estimated radius of 0.99 ± 
0.03 Rj. The black horizontal lines show the corresponding equatorial 
velocity of 36 ± 4 km s"' . Crucially, this equatorial velocity and the as- 
sociated alignment only apply for the period of 3.32 ± 0. 15 hours and a 
radius of 0.99 ± 0.03 Rj. The measured equatorial inclination of 32 ± 4 
degrees is marked with the downward arrow, which is within 10 degrees 
of the orbital inclination angle of 41.8 ± 0.5 degrees (Konopacky et al. 
2012). 



However, an analysis by Dupuy & Liu (20 11), finds that the ec- 
centricity distribution of ultracool dwarf binaries is statistically 
distinct from that of solar type systems, and more consistent 
with the clustered, competitive accretion model of Bate (2009, 
2012). Even if such systems are born misaligned, tidal torquing 
between disks can re-align close systems (see Lubow & Ogilvie 
(2000)). Interaction with a circumbinary disk, as seen in Bate 
(2012), might also align stars with initially random orientations. 

As noted above, magnetic interactions, which might be even 
stronger for fully convective stars, can also alter spin-orbit align- 
ment. The ffcw-dynamo, thought to be responsible for the gener- 
ation of magnetic fields of main sequence stars, no longer func- 
tions for fully convective objects. Magnetic fields many orders 
of magnitude greater in strength than solar-type systems have 
been shown to exist in the low mass star regime (Hallinan et 
al. 2007; Reiners & Basri 2007; Antonova et al. 2008; Berger 
et al. 2009). These fields could have a physical effect on align- 
ment. 2MASS J0746-H20AB, at a separation of only -2.7 AU 
possesses a large-scale ~1.7 kG magnetic field (Antonova et al. 
2008). Indeed, a strongly misaligned system would be more in- 
dicative of dynamical processing. 

Although alignment in 2MASS J0746H-20AB cannot be used 
to distinguish between various formation models, more measure- 
ments of similar systems might elucidate trends with mass. For 



Fig. 6. Equatorial rotational velocity {y-axis, left) vs. inclination 
angle of the orbital plane (x-axis, bottom) for 2MASS J0746+20B. As 
before, the measured v sin i of 2MASS J0746+20B is represented by 
the blue solid curve, where dash-dotted lines represent the errors in each 
measurement. We highlight our estimated radius of 0.96 ± 0.02 Rj (y- 
axis, right), which, together with the radio period of 2.07 ± 0.002 hours 
(Berger et al. 2009), infers an equatorial inclination angle of 36 + 4 
degrees. 



example, one might expect the separation at which systems tran- 
sition from aligned to misaligned to be smaller for lower mass 
systems, particularly if magnetic fields play a role. Comparing 
the orbital properties of stars across the mass spectrum may in- 
dicate where different formation pathways dominate. 



7. Conclusion 

We report on the orbital coplanarity of the L dwarf tight binary 
2MASS J0746+20AB. Recently, high-precision dynamical mass 
and individual rotation velocity measurements were obtained 
(Bouy et al. 2004; Konopacky et al. 2010, 2012), as well a rota- 
tion period for one component (Berger et al. 2009). We present 
the discovery of the rotation period from the primary component 
in this work. From these data, we infer that the binary orbital 
plane is oriented perpendicular to the stellar spin axes to within 
10°. Such alignment has previously been observed in studies 
of solar-type binaries (Hale 1994), and also more massive stars 
(Wheelwright et al. 2011). This work is the first direct evidence 
of spin-orbit alignment in the VLM binary regime. We outline 
the numerous binary formation models that are consistent with 
the observed alignment. Further theoretical work, and a larger 
sample of VLM systems (a campaign that we have already com- 
menced), will place tighter constraints on the most likely forma- 
tion pathways. 

The discovery of the 3.32 + 0.15 hour period implies that the 
radio period of 2.07 + 0.002 hours is associated with the sec- 
ondary star and not the primary star, as was suggested in Berger 
et al. (2009). This in turns refutes the claimed radius of 0.78 
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+ 0.1 Rj for 2MASS J0746+20A. We find that the primary and 
secondary have radii of 0.99 ± 0.03 Rj and 0.96 + 0.02 Rj re- 
spectively. 
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